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High precision membranes are used in many current space applications. This paper 
presents a fully nonlinear membrane theory with forward and inverse analyses of high 
precision membrane structures. The fully nonlinear membrane theory is derived from 
Jaumann strains and stresses, exact coordinate transformations, the concept of local 
relative displacements, and orthogonal virtual rotations. In this theory, energy and New- 
tonian formulations are fully correlated, and every structural term can be interpreted in 
terms of vectors. Fully nonlinear ordinary differential equations (ODEs) governing the 
large static deformations of known axisymmetric membranes under known axisym met- 
ric loading (i.e., forward problems) are presented as first-order ODEs, and a method for 
obtaining numerically exact solutions using the multiple shooting procedure is shown. A 
method for obtaining the undeformed geometry of any axisymmetric membrane with a 
known inflated geometry and a known internal pressure (i.e., inverse problems) is also de- 
rived. Numerical results from forward analysis are verified using results in the literature, 
and results from inverse analysis are verified using known exact solutions and solutions 
from the forward analysis. Results show that the membrane theory and the proposed 
numerical methods for solving nonlinear forward and inverse membrane problems are 
accurate. 

Keywords: Nonlinear membrane theory; numerically exact analysis. 


1. Introduction 

As the cargo space of a launch vehicle is limited, large spacecraft structures 
must be designed to be stowed during launch and deployed once on orbit. 
Moreover, as the launch expenditure of a NASA space mission always consti- 
tutes a significant fraction of the total cost, inexpensive launch vehicles with 
small payload masses are always desirable. Hence, instead of using previous 
electromechanical types of deploying systems, recent efforts of NASA has con- 
centrated on the use of inflatable/deployable structures for space applications . 1 
Inflatable/deployable membranes have been used in space as parabolic antennas, 
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radiators, solar concentrators, sun shields, habitats, radio-frequency structures, 
optical communication systems, radars, lightweight radiometers, telescopes, etc. 
Large balloons are also membrane structures that have been used for many scientific 
missions. 

Materials for membrane structures are usually required to have (1) low 
weight (<0.2 kg/m 2 ), (2) high tensile strength, (3) high tear strength, (4) low 
Helium permeability, (5) appropriate thermal properties, (6) low crease sensi- 
tivity, (7) seam strength greater than or equal to 80% of the parent material 
strength, and (8) no significant losses in tensile or peal strength at low temper- 
ature (e.g., -30°C). 2,3 Materials for current scientific membrane structures include 
reinforced films and coated fabrics. The reinforced films are generally polyethy- 
lene, polyester (e.g., Mylar), and nylon; and the coated fabrics are typically nylon 
or polyester fabric coated with silicone, polyethylene, polyurethane, polyester, or 
nylon. 2 The majority of these materials are being developed by racing yacht sails, 
packaging, geomembrane, parachute, and textile industries. The advantages of 
membrane structures include small stowed volume, light weight, low cost, and good 
thermal and damping properties. 4,5 However, there are difficulties in the design of 
large scientific membranes. 6 

A few factors need to be considered in designing large membrane structures, 
which include construction of smooth membrane surfaces using gores and seam 
tapes, minimum folding requirement for stowage, deployability, accuracy of the de- 
ployed shape, positive principle stresses to prevent wrinkling, large principle stresses 
to provide enough global bending stiffness, small principle stresses to prevent ma- 
terial creeping and change of the deployed shape, and easy manufacturing. These 
are complex, stringent, and somewhat conflicting design requirements. Hence, there 
are unique issues that need to be solved in order to assess and satisfy these design 
requirements. 


1.1. Modeling issues 

Membranes are two-dimensional (2D) analogy of cables, and they have thin thick- 
ness and negligible bending stiffness. Hence, membranes support loads mainly by 
in-plane stretching and, sometimes, shearing. Since stretching is the main load path, 
thickness change due to Poisson’s effect can be significant. For a material with a 
Poisson ratio close to 0.5 (e.g., rubber), the material stiffness matrix is singular if 
the thickness change due to Poisson’s effect is neglected. Moreover, membranes are 
usually designed for undergoing large displacements and rotations but not for small 
or even large strains during regular operations. 2 Large rotations cause geomet- 
ric nonlinearities and make some strain and stress measures nonobjective. 7 Hence, 
choosing an appropriate strain measure that excludes rigid-body displacements and 
rotations becomes the major factor in deriving a fully nonlinear membrane theory. 

Wrinkling is a local buckling phenomenon of 2D structures with negligible 
bending stiffness. A membrane theory without accounting for bending stiffness 
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predicts that wrinkling will happen when the smallest in-plane principal stress de- 
creases from positive to negative. On the other hand, a membrane theory accounting 
for bending stiffness predicts that wrinkling will happen only when the local tan- 
gent stiffness decreases from positive to negative, which usually corresponds to a 
negative but small principal stress. Furthermore, only a membrane theory with 
bending stiffness can be used for post-buckling analyses to predict the bifurcated 
deformation shapes after wrinkling. Although rolling rather than folding is usu- 
ally preferred in packing inflatable structures, 1 folding is necessary in many cases. 
Unfortunately, folding may create permanent folding lines in crease sensitive mem- 
branes, and such folding lines cause surface distortion of the inflated membranes. 8 
To study the influences of folding lines requires the use of a membrane theory that 
accounts for bending stiffness. 

Although the inflated shape of a membrane can be designed by varying the 
distribution of membrane thickness, research shows that temperature change may 
degrade prediction accuracy much more severely than thickness variations. 9 

Hence, a membrane theory for analyzing scientific membranes should account 
for geometric nonlinearities due to large rotations, thickness change due to Poisson’s 
effect, bending stiffness, and thermal influences. 


1.2. Analysis issues 

Analysis of scientific membrane structures is required with an accuracy beyond the 
customary tolerances of conventional structural engineering. For example, a shape 
prediction tolerance of 0.1 mm for a radio-frequency membrane reflector having 
a diameter of a few meters may be taken as the required accuracy of analytical 
solutions. 9,10 

As membrane materials are usually viscoelastic, membrane structures are 
designed to have small principle stresses in order to reduce the possibility of creep. 
However, a membrane with small stresses is vulnerable to wrinkling when the struc- 
ture is subjected to external localized disturbing loads. Because of the nonlineari- 
ties and multiple loading and unloading paths, the wrinkles may stay even after the 
disturbing loads are released. On the other hand, because the magnitudes of inter- 
nal stresses determine the global bending stiffnesses of a membrane, the principle 
stresses should be designed to be positive and as large as creep concerns will allow. 

Although creep is less of a problem for membranes because the internal pres- 
sure gradient can be adjusted and hence creep can be viewed as a change in the 
initial shape. Nonuniform creep due to nonuniform stresses may cause the deployed 
shape to deviate from the required one even after makeup gas is replenished. 
Hence, it is better to design a membrane to have uniform internal stresses after 
deployment. 

The change of deployed shape and redistribution of loads and load carrying 
capabilities of a membrane after wrinkling need to be understood in order to design 
accurate and durable membranes, which requires complex stability and bifurcation 
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analyses. Moreover, because membranes may experience strain values as large as 
14% before burst 2 and membrane materials may be viscoelastic, their material 
properties may depend on strains and/or strain rates and need to be obtained 
experimentally before design analysis. 

Although an initially flat membrane can be inflated to approximate the geometry 
of a parabolic antenna, the inflation pressure required is so high that the makeup 
gas needed to replenish possible loss of gas due to meteoroid damage is also high and 
hence heavy. 11 Hence, it is better to use an initially curved membrane to design an 
accurate deployed shape with a low inflation pressure. However, the design problem 
becomes how to find an unknown undeformed membrane geometry with unknown 
thickness distribution that will result in the required deployed shape with a known 
inflation pressure, which is a challenging nonlinear inverse problem. 

In the past, space systems experienced high failure rates during deployment but 
usually performed well after deployment. 1 With inflatable structures, the stake is 
even higher because their load carrying ability is commensurable with their state of 
deployment. To prevent deployment failure, simulation of the transient dynamics 
of inflatable structures is needed, which requires advanced modeling and analysis 
techniques. 

The above discussions show that designing scientific membranes is challeng- 
ing because of geometric nonlinearities, accuracy requirement, viscoelastic material 
properties, creep, wrinkling, inverse design problems, and transient deployment 
dynamics. However, an accurate nonlinear membrane theory plays the key role in 
solving these problems. 


1.3. Objective 

The objective is to derive a fully nonlinear membrane theory that can be used for 
accurate nonlinear analysis of precision membranes. Moreover, because rotation- 
ally symmetric membranes are popular in scientific applications, nonlinear ODSs 
governing their static deformations under axisymmetric loading will be listed and 
a numerical method for investigating their load-deflection characteristics will be 
derived and numerically verified. Furthermore, a method for computing the un- 
known undeformed membrane geometry that is needed for manufacturing from the 
known/required inflated geometry and the known/designed inflation pressure will 
be derived and numerically verified. 


2. Fully Nonlinear Membrane Theory 

Membranes as well as plates and shells are 2D (two-dimensional) structures. An 
initially curved membrane is a thin shell with negligible bending and twisting stiff- 
nesses. Hence, we first treat a membrane as a doubly-curved shell in deriving a 
fully nonlinear membrane theory for general use, and then we apply assumptions 
to obtain governing equations for simplified types of membranes, especially for thin 
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Fig. 1. The undefonned and deformed geometries of a differential reference surface and the three 
coordinate systems used. 

axisymmetric membranes. The two major tasks in deriving a nonlinear 2D struc- 
tural theory are (1) how to describe the cross section warping, and (2) how to 
describe the reference surface deformation. Different approximations of the cross 
section warping result in different shear deformation theories, and different ap- 
proximations of the reference surface deformation result in different geometrically 
nonlinear theories. As membranes have very thin thickness, the strain energy due to 
cross section warping is negligible. Next we present geometrically exact descriptions 
of the imdeformed and deformed geometries of an initially curved surface. 


2.1. Surface analysis 

2.1.1. Coordinate systems and curvatures 

Two coordinate systems are required in order to describe the deformation of an 
initially curved 2 D structures: one describes the undeformed reference surface and 
the other describes the deformed reference surface. Figure 1 shows an infinitesimal 
area of the reference surface of a 2 D structure before and after deformation. The xyz 
is an orthogonal curvilinear coordinate system with the x and y being curvilinear 
axes on the undeformed reference surface and the 2 being a rectilinear axis, and 
the £7/C is an orthogonal curvilinear coordinate system with the £ and rj being 
curvilinear axes on the deformed reference surface and the C being a rectilinear axis. 
An inertial rectangular coordinate system abc is also used, for reference purpose, 
in the calculation of initial curvatures. We let j 1? j2, and j 3 denote the unit vectors 
along the axes x, y, and z; ii, h, and i3 denote the unit vectors along the axes £, 
77, and Q and i a , i*,, and i c denote the unit vectors along the axes a, 6, and c. 

To obtain initial curvatures, the undeformed position vector P of point A is 
presented as 12 : 


P = Pi (x, y)i« + Pi{x, y)ib + y)i c - 


( 1 ) 
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Taking first- and second-order derivatives of Eq. (1) with respect to x and y yields 


BP 

ji = — -Plxia 4 P2xh 4 , 

dP 

J2 = = Ply^a 4 P2y^b 4 P^y^-c ? 

oy 


(2) 


j3 = jl X j 2 

= ( P2xP^y ~ P3xP2y)ia 4 (i^x^ — PlxP3y)h 4 (-Plx^y “ ^2x^ly)ic , 


and 


^{jl23} 

dx 



' 0 

4 

-4 ' 

= [-K?] {J 123 } 1 

[*?] = 

-4 

0 

~ ^61 



. 4 

^61 

0 

0{jl23} 

dy 


; 

’ 0 

4 

“^62 

= [*§] {jl23> , 


-4 

0 

-k% 



. &62 

4 

0 


where ( ) x = d( )/dx, ( ) y = 3( )/dy, {ji23} = {jij2, j3} r , and [if?] and [K$] are 
initial curvature matrices, /c? and fc 2 are curvatures of bending with respect to the 
axes y and — x, respectively. fe§ x and fc§ 2 are curvatures of twisting with respect to 
the axes — x and y, respectively, k $ is the spiral curvature of the y axis with respect 
to the z axis; and k® is the spiral curvature of the x axis with respect to the z axis. 
The initial curvatures can be derived using Eq. (3). For example, 


k 


0 

1 



~Plxx{P2xPsy ~~ P3xP2y) p2xx{PzxPly PlxPsy) 


~ P3xx{P\xP2y PlxPly) ■ 


( 4 ) 


In Fig. 1, u, u, and w are the displacements of point A with respect to the axes x, 
y, and z, respectively. The axes £ and r) represent the convected configurations of the 
axes x and y, respectively; and 76 (— 761 4 762) is the in-plane shear deformation. 
We let ij and i 2 denote the unit vectors along the axes £ and 17, respectively. We 
note that the axes £ and 77 coincide with the axes £ and fj only if the in-plane shear 
deformation 75 is zero. Letting D = AA / = uji 4 4 wh and using Eq. (3), we 

obtain 

A'B' = -D 4- dxji 4 D 4- ^ dx = dx ji 4- ^ dx 

= [(1 4- u x - vkl 4- wki ) ji 4- (v x 4 uk® 4- wk^) j 2 
+ (w x - uk® - vk^hldx 


(5) 
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A'C' = -D + dyj 2 + D + — dy = dyj 2 + ^-dy 

ay dy 

= [(tty - vk® + w/c6 2 )ji + (1 + Vy + uk° + wk%)} 2 
+ ( Wy - uk% 2 - vk°)j 3 ]dy . 

Hence, the axial strains along the axes £ and r ) , e\ and e 2 , are given by: 

AW-dx 

61 = dx 

= [(1 4- u x - vk° + wki) 2 4- (v x 4- ufc° + wk^) 2 
+ (w x - uki - ^'fc6i) 2 ] 1/2 - 1 7 

Wa-dy 

62 dy 

= [(tty - vk\ + wk% 2 ) 2 + (1 4- Vy + uk° 4- wk 2 ) 2 
+ (w y — uk % 2 - ffca) 2 ] 1 / 2 - 1 • 

The unit vectors along the £ and fj directions are given by: 


(1 4- ei)dx 


= Tnji + Ti2j2 + fi3j3 , 


where 


(1 4- e 2 )dy 


= Tbiji + T22j2 + T23j3 , 


Tii = 


_ 1 4- u x — vk 5 + wk° 
= 1 + ei ’ 

rf, _w x -ukl- vk^ 

-*13 - i 

l + ei 

_ 1 + Vy 4- ufc° 4- wk 2 
1 + e 2 


_ v x + nfcf + wkjji 
1+ei 

u y — vk% + wk% 2 
l + e 2 ’ 


rf, Wy- uk% 2 - Vk% 

J 23 = — 

1 -he 2 


Tl2 = 


?21 = 


It follows from Fig. 1 and Eqs. (8) and (9) that 

76 = 761 + 762 = sin _1 (i| • i^) = sin^TuTbi + T12T22 + T13T23) . (Ha) 

Hence, 75 can be represented in terms of u, v, and w. Moreover, one can use the 
symmetry of Jaumann shear strains to prove that 7 

(1 + ei) sin 7 ei = (1 + e 2 ) sin 7 62 . (lib) 

Then, the unique expressions of 761 and 762 can be determined from Eqs. (11a) and 
(lib). The unit vector normal to the deformed reference surface, i3, is given by: 


= ?3ljl + T32 h 4- T33j 3 , 
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where 


rri rri rri rri rri m rri rjh 

i 12^23 - -L 13-*- 22 ™ ^13^21 “ ^11^23 

= n » J 32 = B 

JLQ Jto 

- ^11^22 ~ ^12^21 
i?0 


( 13 ) 


Ro = [(^ 12^23 “ T13T22) 2 + (T13T21 — T11T23) 2 + (fiif 22 “ r^Tbl) 2 ] 1 ^ 2 


= | COS 76 1 . 


Using Fig. 1 and Eqs. (8), (9), and (12), we obtain the following transformation 
which relates the undeformed coordinate system xyz to the deformed coordinate 
system 


where 



'tn 

T 12 

t 13 ' 

{ii23} = [T]{ji 23 }, [T\ = [T] 

t 2 1 

T22 

tj 3 


T z 1 

T32 

7 33. 


{il23> = [r]{i|2 3 } 


[T] 


COS 761 

sin 762 

0 


sin 761 0 

COS 762 0 

0 1 


1 

cos 76 


COS 762 

- sin 762 

0 


- sin 7 ei 

COS 761 

0 


0 ' 
0 

cos 76 . 


(14a) 


(14b) 


{»123} = {ii, 12 , h} T i and {ij^ 3 } = {ij,i 2 ,i 3 } T . Moreover, using Eqs. (14a) and (3) 
and the following identities 


di j . d\j . dij . di k . 

- ■ >* = 


we obtain 


^{ 1123 } 

dx 


dy 


= [^l]{ll23} , 


dy 


[*i] = 


0 k 5 — k\ 

— 0 — &61 

L fei fc 61 0 J 


= ^[T} t + [T]{K° 1 }[T\ t , 


(15) 


(16) 
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d{im} 

dy 


[K 2 } = 


= [# 2 ] {> 123 } , 
0 * 4 
— *4 0 

*62 *2 


*62 

*2 

0 


= ^B [ T}T +[ t}[K^\T} T - 


(17) 


where [K\] and [K 2 ] are called deformed curvature matrices in this work and 


01 

kl = * * 3 = ~T\rnxTz m — + 322^1 + I23&5 , 

*2 = ~ ■ i 3 = ~T 2my T 3m + T n k° - J 12 *g 2 - r 13 *? , 
*61 = ~ ■ i3 = -T 2 mxT 3m + T n k §! - r 12 *? - 7i3* 5 ° , 

*62 = ' * 3 = -T\myT 3m ~ 721*2 + 7 22 *§ 2 + ^23*4 i 

*5 = ^ • J 2 = r lmx r 2m - t 31 * §1 + r 32 *? + Tasfcg , 

*4 s • ii = -T 2my T lm - T 31 k% + T 32 *g 2 + T 33 * 4 ° . 


(18) 


Here, ii = i 2 x i 3 , i 2 = i 3 x ii, 13 = ii x i 2 , and Eq. (14a) are used. Hereafter, 
unless otherwise stated, repeated subscript indices imply summations. We note 
that the kj in Eq. (18) are functions needed for exact description of £ and 77, but 
they are not real curvatures because the deformed dx (dy) is not along the ii (i 2 ) 
direction because 76i ^ 0(762 ^ 0). If 761 = 762 = 0, the curvatures are normalized 
(but not real) curvatures because the differentiations in Eq. (18) are taken with 
respect to the undeformed lengths dx and dy, instead of the deformed lengths 
(1 + ei )dx and (1 + e 2 )dy. Only if 761 = 762 = ei = e 2 = 0, kj represent real 
curvatures. 


2.1.2. Variations of global strains 

To derive the governing equations using the extended Hamilton’s principle, we 
need to have the variations of the global strains ei, e 2 , 7 6 i, 762> and ki. Taking 
variations of Eqs. (6) and (7) and using Eq. (10), we obtain the variations of e\ and 
e 2 as: 


Se 1 = TuStn + Ti 2 5ti 2 + Ti3<5ti3 , (19) 


£e 2 = T2l<ft 2 l + ^22^22 + ^23^23 > 


( 20 ) 
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where 


<ftii = ^(1 + %- vk 5 4- wk ^ ) = 5u x - k%5v 4 , 

toi2 = 5{v x 4 U&5 4 icfe^) = 5i/ x + k®6u 4 , 

toi3 = tifc? - vfeei) = Su) x - kiSu - k% x Sv , 

to 2 i = 5(w y - 4 tufc&a) = Su y - 4 fcg 2 5ie , 

to 2 2 = 5(1 + v y 4- ix/c$ 4- wfc 2 ) = $ v y + 4 fc 2 5w , 

5^23 = 5(tu y — uk $ 2 — v& 2 ) = 5w y — k^ 2 5u ~ k 2 5v . 

Moreover, taking variations of Eqs. (8) and (9) and using Eq. (10) 

yields 

5ij = — — - — (jiton 4 j 2 5ti2 4 js5ti3 — i^ei) , (22a) 

i 4 ei 


5i 2 — ~ (ji5t 2 i 4 j 2 5t 22 4 J35t 2 3 — i 2 5e 2 ) . (22b) 

Taking the variation of sin 7 6 = ij • (see Eq. (11a)) and using Eqs. (22), (8), (9), 
(19), and (20), we obtain 

_ (T 2 i - sin 76Tn)5tii 4 (T 22 — sin 76Ti 2 )5tj 2 4 (X23 — sin 76^13)5^13 

cos 76(1 4 ei) 


{Tn — sin 76T 2 i)5t 2 i 4 (Ti 2 — sin 76T 22 )5t 22 4 (T13 ~~ sin 76^23)5^23 

cos 76 (1 4 e 2 ) 


(23) 


Taking the variation of Eq. (lib) and using the fact that 576 = <$761 4 5762, we 
obtain 


5761 = 


(1 4 e 2 ) cos 762576 - sin 76i5ei 4 sin 76 2 5e 2 
(1 4 ei) cos 76i 4 (1 4 e 2 ) cos 7 62 


(24a) 


(14 ei) cos 761 576 4 sin 7ei5ei - sin 7 62 5e 2 

0 762 = 7- 7 7- r , (24b) 

(1 4 ei) cos 761 4(14 e 2 ) cos 7 62 

where 576 is shown in Eq. (23). 

To obtain the variations of curvatures (i.e., Skj), the concept of orthogonal 
virtual rotations is needed. Because the variations of the unit vectors are due to 
the virtual rotations of the observed membrane element, we have 


0 503 —50 2 

—503 0 50i , 

50 2 -50i 0 


<5{»123} = [<50]{il23} , [<5#1 = 


(25) 



Fully Nonlinear Modeling and Analysis of Precision Membranes 29 


where 60\, S0 2 , and 503 are the virtual rigid-body rotations of the observed mem- 
brane element with respect to the axes £, r/, and (, respectively. We note that 80i are 
differential rotations and hence they are vector quantities. Moreover, 56 1 are along 
three perpendicular directions and hence they are mutually independent. Taking 
the variations of curvatures defined in Eq. (18) and using Eqs. (16), (17), and (25), 
one can obtain that 12 


' - 5*61 ' 

< Ski ► 
k 5k 5 i 

'-Sk 2 'j 
• Ske 2 > 

k Sk 4 J 


d_ 

dx 


d_ 

dy 


'SO,' 

r 

so 2 

' - [Ki] < 

JO 3j 


'<501 j 

r 

<5^2 

\-[K2}{ 

. <503 J 

{ 


69 2 




f *1 


> 


(26a) 


(26b) 


Using Eqs. (25), (14b), and (22) and the fact that i 3 - = i 3 - — 0, we obtain that 


66i = Si 2 • i 3 = • is - -is 


COS 76 


cos 76 


cos 761 


cos 76 (1 4 e 2 ) 


(T 31 St 21 4- T^ 2 St 22 4 Tz^Stys) 


\ C^n 4 T 32 5fi2 4 T&6t 13 ) , (27) 

cos 7 6 (1 4 ei) 


S6 2 — — 5ii • i 3 — 


sin 76i 

COS 76 


5i2‘i 3 - 


COS 762 
COS 76 


5ii * i 3 


_ ^ 61 — (r 3 i5t 2 i 4 T 32 5^22 4 TssSt 2 3) 

cos 76(1 4 e 2 ) 

— _ (r 3 i5t u 4 T 32 5fi2 4 T 33 5ti 3 ) . (28) 

cos 76(1 4 ei) 

Using Eqs. (25), (14b), and (22), we obtain that 
SO 3 = \ (fli • i 2 - <5i 2 • ii) = l (*76 2 - *m) + x-4- (Si i • i 2 - Sii ■ ij) 

Z Z Z COS 76 

= 5 (<5762 - <5761 ) 

(jT 2 i — sin 26fn)Stn + (J 22 — sin 76Ti2)<5ti 2 + (!T 2 3 — sin 76ri3)<5ti3 

2 cos 76(1 + ei) 

_ (Tii — sin 76^21)^21 + (fui ~ sin 76r 22 )^22 + (T~i3 - sin 76^23) <5^23 

2 cos 76(1 + e 2 ) 


( 29 ) 
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Hence, be i, be 2, £75, 6761* <^762> J0i, J02> and b9^ can be represented in terms of 
Jtt, bv, bw , Ju x , Ju x , <5u> x , Ju y , bv y , and Jit; y . Moreover, the variations of curvatures 
(Jfcj can be represented in terms of Ju, du, <5u?, 5u x , bv x , Ju; x , <Ju y , Ju y , <5u; y , <5u xx , 
*^xxi ^^xyi ^^xy> and blVxy • 

2.2. Jaumann strains and stresses 

Although membranes have thin thickness, it is necessary to account for bend- 
ing stiffness in order to study wrinkling problems and shape control of inflatable 
membranes. For the modeling of highly flexible structures, Jaumann strains and 
stresses are appropriate and convenient measures because they are objective, work- 
conjugate, and fully nonlinear. 7 Moreover, because Jaumann strains are corotated 
engineering strains, material properties obtained from experiments using small en- 
gineering strains and stresses can be directly applied. 

Next we derive the fully nonlinear expressions of Jaumann strains by using the 
polar decomposition theory and the concept of local relative displacements. 


2.2.1. Polar decomposition 

For an arbitrary point Q on a differential membrane element that has dimensions 
dx x dy x h (h = thickness) and crosses the reference surface at point A shown in 
Fig. 1, because transverse shear warping is assumed to be negligible, the undeformed 
position vector P q and the deformed position vector P q are given by: 

Pq = xh + yh + zj 3 , (30a) 


P Q = (x 4- u)ji + (y 4- t))j2 + u>j3 + 2(1 + 63)13 (30b) 

where e% is the extensional strain along the thickness direction and is assumed to be 
independent of 2. Because x and y are local curvilinear axes starting from point A, 

x = y = 0. (31) 


Taking derivatives of Eq. (30a) and using Eqs. (3) and (31), we derive a gradient 
tensor [Y] as: 


F] = lYmn} M 


dP Q ' 


r* dx n \ 



1 + zki 

zk 61 

0 


zk$ 2 0 " 

1 4- zk% 0 , 

0 1 


(32) 


where x n (n = 1,2,3) denote x, y, and 2, respectively. Similarly, taking derivatives 
of Eq. (30b) and using Eqs. (3), (16), (17), (14a), (31), and (10), we obtain the 
deformation gradient tensor [Y j as: 
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[Y] — [Vmn] = jm * j 

6 iXn 4 - z(kiTn + £61X21) + ze^xTzi S2T21 4 - £(fce 2 Xn 4 - £2X21) 4 - ze 3 yX 3 i £3X31 

61X12 4 - z(k\T\2 4 - £61X22) 4 - *e 3 x X 32 62X22 4 - £(£62X12 4 * £2X22) 4 - zesyT^ 63X32 

61X13 4 - £(£1X13 4 - £61X23) 4 - ze^xT^ 62X23 4 - £(£52X13 4- £2X23) 4 - ze 3 yX 3 3 €3X33 

(33) 

where 2 = (1 4- 63 ) 2 , ei = 1 4- ei, £2 = 1 -f e 2 , and e 3 = 1 + e 3 . It can be shown by 
using the polar decomposition theory that 13 

[Y] = [T} t {U], (34) 

where [U] is the so-called right stretch tensor and should be a symmetric matrix, 
and [T] accounts for the corotation from the system xyz to the system £ 77 ^ as shown 
by Eq. (14a). Jaumann strain tensor \B] is defined to be 

[B] = [[/]-[/], (35) 

where [/] is a 3 x 3 identity matrix. It follows from Eqs. (34), (15), (33), and (14a) 
and the orthogonality of i, that 

' ex cos 76 i + zki e 2 sin 7 62 + zk 62 0 

[17] = [T\ [Y] = ex sin 7 6 x + zlc^i e 2 cos 7 62 + zk 2 0 

_ e 3 

To determine principal strain directions of [ B ], one needs to determine the eigen- 
values A of [U] by solving 

|[C/]-A[/]| = 0. (37) 

Equation (36) indicates that the eigenvalues A* will vary with 2 . However, because 
the membrane thickness h is usually small, one can choose the reference point 
(z = 0) for computing principal strain directions. The eigenvalues of [U ] z = 0 can be 
obtained to be 

Ai, A 2 — 

A 3 = 1 4- e 3 . (38) 

The corresponding eigenvectors are: 

{Ci} = {Cii,C 2 i,0} 7 \ {C 2 } = {— C 21 ,Cxx,0} t , {C 3 } = {0,0, 1} T , (39) 

where Cu and C 2 i can be obtained from Eq. (36) (with 2 = 0) as functions of ei, 

C 2 , 761 , and 762 * If 761 — 762 = 0, eigenvalues of [U] z = 0 are 1 + ei, l + e 2 , and l + e 3 , 
and C\\ = 1 and C 2 i = 0. In other words, the axes £, 77 , and Q are the principal 
strain axes. Pai and Palazotto 13 showed that, if {(7,} are normalized such that 

[Cf[C] = [/], (40a) 


ex cos 761 + e 2 cos 76 2 ± \J(e\ cos 7 6 i + e 2 cos 762 ) 2 - 4eje 2 cos 7 6 

2 
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then 

" Ai 0 0 " 

{ii23} = [<?] T {ii23}, [ U] = [C] 0 A 2 0 [C] T , (40b) 

.0 0 As. 

where [C\ = [{Ci}{C 2 }{Cz}}, {ij^} = {h h h} T > an( ^ h are base vectors along the 
principal strain directions. 

Equations (lib) and (36) show that [U] z =o is symmetric. For strains at z / 0, 
[U] is symmetric only if k^i = k $ 2 and e3 X = e% y = 0. ke 1 = ^2 is true if the 
deformed reference surface is smooth. However, e3 X and e$ y may not be zero and 
approximations due to computation may make ke 1 7^ &62? and hence the [U] in 
Eq. (36) is usually asymmetric. The transformation matrix [T] in Eq. (14a) is 
derived only for points on the reference surface (i.e., z = 0), but it is used in 
Eq. (36) for all points. Hence, the [U] in Eq. (36) is not really the right stretch 
tensor for points at z ^ 0 because of the nonuniform distribution of e3 and the 
local rotations due to twistings (i.e., £fc 61 , £^62). Consequently, instead of using 
Eq. (35), Jaumann strains I? mn are modified as: 

[B] = \[{U} + [U} t -[Y\-[Y] t ] . (41) 

Substituting Eqs. (32) and (36) into Eq. (41) yields the Jaumann strains of a doubly- 
curved membrane as: 

£11 = (1 + ei) cos 761 - 1 4- zki - zk J , B 22 = (1 4- e 2 ) cos 762 - 1 4- zk 2 - zk \ , 

Bn = \ [(1 + ei) sin 7ei + (1 + e 2 ) sin 7 62 + zk 6 - zk %] , (42) 

i 1 

-D33 = C3 , B 13 = - e 3x Z , B 23 = ^ e 3 y z , 

where k$ = &6i +&62 and = k ^ 4- ^62- As membranes usually have thin thickness, 
the ^ is small. Hence, we replace £ with z, neglect e^ x z and e$ y z, and rewrite Eq. (42) 
as: 

(B} = [S]W, B 3 3 = e 3 , (43) 

where 

{ B n ] ["I 0 0 2 0 0" 

B 22 ► , [S]= 0 1 0 0 2 0 , 

2 Bi 2 j _0 0 1 0 0 2_ (44) 

W = {(1 + ei) cos 761 - 1, (1 + e 2 ) cos 702 - 1, (1 + ei) sin 7 61 
+ (1 + e 2 ) sin 7 62 , k\ - k\ , k 2 - fc 2 , k 6 - k%} T . 
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2.2.2. Concept of local relative displacements 

Next we show that Jaumann strains can be easily derived using the concept of local 
relative displacements. 14 The local displacement vector u, of an arbitrary point on 
a differential membrane element, has the form 

U = Ufcifc , (45) 

where 


ui(x,y,z,t) = Ui(x,y,t) + z[9 2 (x,y,t) - 0 2 (x,y)], 

u 2 {x,y,z,t) = u°(x,y,t) - z[0i(x,y,t) - 9^{x,y)} , (46) 

u 3 (x, y, z, t) = u£(x, y , t ) + ze 3 . 

Here, u® are the displacements (with respect to the system £r/£) of the reference 
point A! shown in Fig. 1, 9\ and 62 are the rotation angles of the observed membrane 
element with respect to the £ and rj axes, respectively, and 0° and axe 
corresponding initial rotation angles. Because the £t/( is a local coordinate system 
attached to the observed membrane element, and the £ — 77 plane is tangent to the 
deformed reference surface, we have 


It follows from Fig. 1 and Eqs. (3), (16), and (17) that 
du^ _ . du® 


duo dv£ 


dy 


-0. 


(47) 


(1 -f ei)dx 


= sm 7 6 i 


du 1 
dx 


= (1 + ei) cos 761 - 1 , 


= sm 762 , 

(1 4 e 2 )dy 

= (14 e 2 ) cos 762 - 1 , 


du® 

dy 


M = . 5 , = *0 

dx dx J3 fcl ’ 


0 M-_fc 0 M- fc o 

dy ~ 2 ’ dx ~ 61 ’ dy ~ 62 ’ 


(48) 


ai. 


•13 


a^i 


ae 2 


= &62 • 


do 2 aii . a#i <7(7i 

- 5 — = -- 5 — • 13 = «1 , - 5 — = — - 5 — = -«61 1 o 

ax ax ay ax dy 

Taking the derivatives of Eqs. (45) and (46) and using Eqs. (47) and (48), we obtain 

au 

dx 

= [(1 + ei) cos 761 - 1 + zki - zfcj]ii + [(1 + ei) sin 7 6 i + zk^i - zk% i]i 2 
+ ze 3x i 3 , (49) 


dui du 2 du 3 dii di 2 di 3 

ax ax dx dx dx dx 


^ = [(1 + e 2 ) sin 76 2 + zk % 2 - zfc£ 2 ]ii + K 1 + e 2) cos 762 - 1 + zk 2 - zk%]i 2 
dy 


+ ze 3 y i 3 , 


(50) 


au 

dz 


= €313 . 


(51) 
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In terms of the local relative displacements shown in Eqs. (45) and (46), Jaumann 
strains can be simply defined as 14 : 

if du , 9u . \ 

Bran ^ I Qj, * "h * 1?Tl I 1 ft — 1, 2, 3 . (52) 

Substituting Eqs. (49)~(51) into Eq. (52) yields the same Jaumann strains shown 
in Eq. (42). 


2.2.3. Constitutive equation 


For isotropic materials the constitutive equation that relates Jaumann stresses J mn 
to Jaumann strains £? mn is given by: 


' J 11' 



"1-1/ 

V 

V 

0 


' B n ' 

J 22 

i 

E 

V 

1-1/ 

V 

0 

i 

B 22 

J 33 

r 

(1 + i/)(l — 2i>) 

V 

V 

\-v 

0 

B 33 

. J 12 > 



0 

0 

0 

l-2u 

2 J 


, 2S 12 ; 


where E is Young’s modulus and v is Poisson’s ratio. We note that, for a material 
with v = 0.5 (e.g., rubber), numerical singularity exists in the constitutive equation. 
However, because the thickness of a membrane is usually thin and hence £? 33 (= e 3 ) 
is mainly due to Poisson’s effect, one can assume that J 33 = 0 to obtain that 

^3 = t Bn H B 22 * (54a) 

Consequently, we have 



771 

"1 

V 

0 ■ 


' B n ' 

b 

E 

k — 

V 

1 

0 

i 

B 22 

{ 2-Bi2 > 

1-^2 

0 

0 

1 — V 

2 J 


(54b) 


Membrane materials can be anisotropic. For example, Mylar has been reported 
to be orthotropic. 11 If the membrane is made of layers of orthotropic materials, 
one can obtain the transformed stiffness matrix [QW] for the ith lamina from its 
principal stiffness matrix [QW] and its ply angle (measured with respect to the axis 
x) by using tensor transformation and assuming J33 — 0. Then, the constitutive 
equation for the zth lamina is given by: 

{J} = [Q«]{B}, (55a) 

where 



f jd) } 


ros 

Ql2 

OS' 

{J}=< 

t(*) 

J 22 


OS 

Q22 

Q 26 


rW 

l J 12 J 


Los 

Q26 

@66 m 


(55b) 
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Moreover, stress resultants are defined as: 

j(0 

J n 
J 22 

f 12 

h zJ$ 

Z ^22 

Ml 

\ [0 (i) ][5] I f \ [<3«] z[QV] 1 

^[g (i) ][S]J Jz [z[Q®] z 2 [Q0)]J 

where Eqs. (43) and (44) were used and [Q] is a 6 x 6 matrix. Here Mi, M2, and 
Mq represent the moment intensities acting on the edges of a differential membrane 
element, N\ and iV 2 are the in-plane extension force intensities, and N$ is the 
in-plane shear force intensity. 

2.3. Governing equations 

To derive governing equations we use the extended Hamilton’s principle, which 
states that 15 

(6T-SU + 6W nc )dt = 0, (57) 

where t is the time, ST is the variation of kinetic energy, (HI is the variation of 
elastic energy, and SW nc is the variation of nonconservative work. 

Deformations of membranes are mainly caused by rigid-body displacements and 
rotations and in-plane stretching, and the in-plane shear strain 76 (= 761 4- 762) is 
usually small. Hence, the influence of 76 on the deformed geometry can be neglected 
without significant loss of accuracy in accounting for geometric nonlinearities. Sub- 
stituting 761 = 762 —76—0 into Eqs. (14a) and (14b) yields 

[r] = [1] ? Tii = Tu , T2i — T2i 7 ii = i| , i 2 = *2 * (58) 

Substituting 7 6 i = 762 = 76 = 0 and Eq. (58) into Eqs. (19), (20), (27), and (28) 
yields 

Se 1 = T\\S(u x - vk$ 4- wk®) + TnS(v x 4- uk 5 + 

4- TizS(w x - uk ? - vk % x ) , (59a) 






Se 2 = T 2 \S{u y - vk\ 4- wk% 2 ) 4- T 2 2 S(v y 4 uk% 4- wk 2 ) 
4- T 23 6(w y - uk% 2 -vk%), 


(59b) 
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-(1 + e 2 )60i + T 31 (8u y — k°Sv + k® 2 Sw) + T 32 (5v y + k®8u + k 2 Sw) 

+ T 33 (8w y - k% 2 Su - k 2 Sv) = 0 , (60a) 


(1 + ei)<S02 + T 3 i(5u x — k®5v + k 2 8w) + T 3 2 (Sv x + k°8u + k^fiw) 

+ T 33 (8w x - k°5u - k^Sv) = 0 . (60b) 

Using Eq. (58) and taking the variations of ij (= ij) and i 2 (= ij) in Eqs. (8) and 
(9), we obtain 

£ii = + j 2 <5Ti2 +j 3 STi 3 

= [(8u x — k 3 Sv + kx8w)ji + (8v x + k®5u + kf-x8w)j 2 4- (Sw x — k^Su — /cgi<5i;)j 3 ] 


5e i . 


1 ■+■ C\ 1 + ej 


(61a) 


^*2 — jl ^21 + J 2^22 + j 3^23 

= [( Su y — k® Sv 4- kQ 2 fiw)ii 4- {$Vy 4- k®6u 4- k^Sw)^ 4- {Sw y — k^ 2 Su — 
1 5e 2 . 


I +€2 l + e,‘ 2 ' < 6lb > 

Using Eqs. (61a) and (61b) and considering 751 and 7^2 being small, we obtain 
that 

(J761 = Sii • i2 

= rZ — P 21 (Su x - k®5v 4- k®6w) 4- T 2 2{Sv x + k®6u 4- 

J- 1 Cj 

+ T2 3 (8w x - k°Su - k^Sv)] , (62a) 


<*762 = <^i 2 • ii 

= — [Tii(8u y — k®8v + k^ 2 8w) 4- T l2 (5v y + k®5u + k 2 Sw) 

+ Ti 3 (5w y - k$ 2 5u - fc^v)] . (62b) 

Moreover, when 76, 761 and 76 2 are small, the Jaumann strains in Eq. (44) reduce 
to 

B n = ex + z(kx - kf ) , B 22 = e 2 + z(k 2 - k 2 ) , 

2B12 = 761 4- 76 2 + z(ke — k %) . 

It follows from Eqs. (30a) and (30b) that the displacement vector D of an 
arbitrary point on the differential membrane element is: 

D = p « ~ Pq = «ji + uj 2 + wj 3 + zi 3 - zj 3 . 


(64) 
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Taking the variation of Eq. (64) and using Eq. (25) yields 

ST> = ji Su 4 j 2 <Sn 4- j 3 5iu + z(\\502 — 12SO1) . (65a) 

Taking time derivatives of D and using Eq. (14a) yields 

D = uji + i>j2 + wj 3 4 zifs-ai 4 1^2 4 T33J3) . (65b) 


Using Eqs. (65a), (65b), and (14a), we obtain the variation of kinetic energy as: 


ST 


-w 


pD * S&dxdydz 


= — J (A u 5u + A v 6v + A w 5w + Aq^Oi + Ag 2 602)dxdy , 


( 66 ) 


where A is the area of the undeformed reference surface, p is the mass density per 
unit volume, and 


A u = mu 4 Jif 3 i , A v = mv 4 /1T32 , A w — mw 4 /1T33 , 

Aq x = —I\(uT2i 4 t)T22 4 WT22) — ^2(^31 -I21 4 T32T22 4 T33T23) , 

(67) 

= h(uTn 4 i)Ti2 4 it>7i3) 4 hi^iTn 4 132T12 4 T33T13) , 
{m,Ji,/ 2 } = / 2 p{l,£,£ 2 }dz, 

where m represents the mass per unit area and z — (1 4 e 3 )z. If p is symmetric 
with respect to the middle reference surface, I\ = 0. Moreover, we also assume 
that the membrane thickness h is small and hence the rotary inertia density I2 is 
negligible. Substituting 1\ = I2 — 0 into Eq. (66) yields 


ST 


— I ( muSu - 
Ja 


mvSv 4 mw5w)dxdy . 


(68) 


If the membrane is only subjected to a normal pressure p along the £ axis, the 
nonconservative work is given by: 

<SWnc = J p(l + ei)dx(l + e 2 )dy\ 3 ■ <5(uji + uj 2 + wj 3 ) 

= f p{T 3 i6u + T 32 Sv + T 33 Sw)(l +ei)(l + e 2 )dxdy 

Ja 

= / (pi<5u 4 4 ps5w)dxdy , (69) 

JA 


where Eq. (14a) is used and the actual deformed area (14 ei)<£r(l 4 e 2 )dy is used 
to calculate the actual load. Moreover, 


pi ee p( 1 4 ei)(l 4 e 2 )T 3 i , P2 = p(l 4 ei)(l 4 e 2 )r 32 , 
p 3 =p(l4ei)(l4e 2 )T 33 . 


( 70 ) 
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Because the elastic energy II is due to relative displacements among material 
particles, SU is the virtual work done by the internal forces f* through the virtual 
relative displacement Su . Hence, we have 


SU 


= [ fi-inS ( T r— dx 1 • i n ) + f 2 • i„£ ( dx 2 • i„ 

Jv° l \dxi ) \dx 2 


( 9u 

+ k-i n S ( — dx 3 *i n 


(71a) 


where dx \ = dx , dx 2 = dy , dx 3 = dz, dF° = dxdydz , and f* is the force vector 
acting on the convected area of dx m dx n (i ^ m ^ n). Because Jaumann stresses are 
defined as 13 : 


dn 


fi-ii 

dx2dx 3 


J22 = 


f*2 • 12 

dxidxs ’ 


dl2 


fi • 12 
dx2dx 3 


f*2 • ii 
dxidx3 


= ^21 , 


(71b) 


substituting Eqs. (71b) and (52) into Eq. (71a) and using B i3 = B 23 = J 33 = 0 
yields 


SU = J J {JuSBii -{- J 22 SB 22 ~h 2di2^^i2)dxdydz . (72) 

Substituting Eq. (63) into Eq. (72), using Eqs. (26a) and (26b), and integrating by 
parts we obtain that 


dll = J J [Jn(Sei + zSk\) + ^ 22(^2 -h z5k2) + di2(<^76l d - ^62 d - zSk§)]dxdydz 
= J (Nide 1 + N 2 Se 2 + N 6 Sj 6 i + N 6 S^62 + M\5k\ + M 2 8k 2 + M 6 8k e )dxdy 
= J (NiSei + N 2 6e 2 + N^ei + + 0i<5#i + Q 2 60 2 + Q^SO^dxdy 

+ f {-M 2 86 1 + dx + [ [-M 6 50 i + M 1 86 2 } x x z^dy , (73) 

J x «/ y 


where X and Y denote the integration limits of the curvilinear coordinates x and 
y , respectively, and 


©i = Mq x + M'2 y + M\ks + M^k^ , 

©2 = - Mi x “ M6y + M2 / C4 + M§k§ , 
© 3 = Mi &6 i — M2&62 + MgA ;2 — Mg/ci . 


(74a) 
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Moreover, it follows from Eqs. (63), (44) and (56) that 


’N x ' 


ei 

n 2 


e 2 

N, 6 


76 

Mi 

- A:® 

m 2 





M-kl, 


(74b) 


Substituting Eqs. (59a), (59b), (62a), and (62b) into Eq. (73) and adding the 
integrals of Q 2 x Eq. (60a) and Q 1 x Eq. (60b), we obtain 

sn = j A [T}{6uvw} x + 

+ |^^,iV 2 ,Q 2 J [T]{6uvw} y +{j^,N 2 ,Q 2 \ {T}[K»} t {6uvw} 

4- [©i — Q 2 (l + 62 )]^i + [0 2 + Qi(l + e x )]S9 2 + ©3J03^dxdj/ 

+ f [-M 2 <50, + M 6 S9 2 } y v zldx + [ [~M g S6 i + M } 50 2 }*z*dy , (75) 

Jx Jy 

where {uvw} = {u v w} T . It is shown later that the introduced Q i and Q 2 are 
transverse shear force intensities. 

Substituting Eqs. (68), (69), (75), and (74a) into Eq. (57), integrating by parts, 
and using Eqs. (60a) and (60b) for the boundary terms — M$86\ and M§802 in 
Eq. (75), we obtain 

L { {~ F dx [T] - + d[F 0y [T] - {F 0 } t [T}{K°} t ^ {6 uvw} 

+ ({i? F } r - {/ F } T ){<W} + ({Q q } t + {Q 0 } T ){50} 


+ 


f d{M a y 
V dx 


+ (m q } t [a: 1 ] + 


d[Mgp_ 

dy 


{M 0 } t [K 2 }^ {69}}dxdy 


~ ( \{F 0 } t {T]{6uvw} + K^ 1 ) + (Tssfc? - T 32 fc°) 1 Su 

Jx I LVl + ei/ x 1 + ei J 

IY M,T 32 \ Me 


[yT^) x + T^-J T ^ +n ' k °\ 


8v 


+ [ (¥^) - TT~ ( r »*6i + T »l*l)l 5w ~ M 2^1 

[ \ 1 4- ei / x 1-fei 




dx 


(76) 


y=0 
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-l 


( MqT^\\ Mq 


Vl + e 2 y u + 1 + 


{F a } T [T}{Suvw} + 

+ \i.TrS) t + ^ iT33k ° +Tnk ° ) 


62 


(T 33 fc° 2 - T 32 fc 4 °) 


Su 


Sv 


4- 


+ 


Mq + Me 


1 + e\ 1 4- e 2 


(5 k; 4- Mi 602 

(T3i(5u 4- Tz2&v 4- T33(5k;) 


x=X 


dx 


x=0 

1 (*,y)=(0,0),(X,V) 


= 0 , 


J (i,j/)=(x,o),(o,y) 


where 


(77) 


{0} = {0 1 ,0 2 ,0 3 } 7 \ 

{/ f } = {mii, mi), rmi;} T , {R F } = {pi,p 2 ,p 3 } T , 

{M q } = {-Me, M x , 0} r , {M^} - {-M 2 , M 6 , 0} T , 

{Qa} = {0, -Qi(l + ei), A^} t j {Q/?} = {<22(1 4- e 2 ), 0, —N G } T 


Setting the coefficients of Su , 5k, 5k;, 50i, 50 2 , and 503 in Eq. (76) equal to zero, 
we obtain the equations of motion as: 


-t (PT{F Q }) - MTf{F Q } + £ m T {F 0 }) - [K o 2 ][T] T {F 0 } + {Ff} = {/ F } , 


d 


dx 


dy 


(78) 


9{Ma} + [Fi] t {M q } + + [K 2 ) t {M 0 ) + {Q q } + {Q 0 } = {0} . (79) 


dx 


dy 


The boundary conditions are to specify: 
Along x = 0, X: 


Su = 0 or F a • ji 4- 

Sv = 0 or F a • j 2 4- 

Sw = 0 or F a • j 3 4- 

50 2 =0 or Mi . 


MqTsi 


1 4- e 2 
MqTz2 


4- 




1 4- e 2 


+ 


J y 


1 + e 2 


J y 


(r 33 fcg 2 - r 32 fc$)M 6 
1 + e 2 

(r 33 fcg + T 31 fcO)M 6 
l + e 2 

(r 32 fc° + r 31 fc° 2 )M 6 

1 + e 2 
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Along y = 0, Y: 

Su = 0 or F# • ji + 

5v = 0 or • j 2 + 

5w = 0 or F^ • j 3 + 
50i =0 or M 2 . 


[ M 6 r 3 i i (jyj - Tyjm 

[l + e 1 \ x 1 + e! 

r M 6 r 32 ] (Taafcg, + r 31 *g)M 6 

[l + e 1 \ x l + ei 

M 6 r 33 l (T 31 fc? + r 32 A&)M 6 
_l + e 1 J I l + ei 


At (z, y ) = (0, 0), (0, r), (*, 0), (X, Y): 


Su = 0 or 
5v = 0 or 
Sw = 0 or 


m 6 

1 

M 6 1 

1 4- ei 

~r- 

1 4-e 2 

m 6 

i 

Me ' 

l + ei 

T 

1 4-e 2 

M 6 

1 

Me ' 

1 4- ej 

I 

1 4- e 2 


7k, 

732, 


733, 


where 


F a 


Mix 4- 


Ne 

1 +ei 


*2 4- Ql>3 , 


F/3 = 


Ne 

1 + e 2 


ii 4- iV 2 i 2 + Q 2 i 3 . 


We also define 


M„ = -M 6 ii + Mii 2 , 

M/j — — M 2 ij 4- M 6 i 2 , 

I F = miijj 4- muj 2 4- mt/>j 3 , 

Rf = Pljl + P% }2 + P3j3 • 


(80) 


(81a) 


(81b) 


Using Eqs. (16) and (17) and the identity [T’] r = [T] x , we rewrite Eq. (78) as: 


[T] T + [Ki ] T {F a } + ^1 + |*m } ) + { Rf } = {If} ■ (82) 

Using Eqs. (81a), (81b), (16), and (17), we put Eqs. (82) and (79) in the following 
vector forms: 


dFq dF 0 

dx dy 


4- Rf = If , 


(83) 


^ + ^ + (l+ei)ix xF Q + (l + e 2 )i 2 xF^ = 0. (84) 

dx ay 

Equations (83) and (84) state the balance of forces and moments on a differential 
membrane element, which can be directly obtained using the stress resultants in 
Eq. (74b) and a vector approach with Newton’s second law. Although we started 
with the energy formulation (see Eq. (57)), the results are fully correlated with 



42 P. F. Pai & L. G. Young 


the Newtonian formulation and every structural term is interpretable in terms of 
vectors. This shows that the equations obtained are intrinsic. 

One can see from Eq. (81a) that Q i and Q 2 are along the ^-direction and hence 
represent the transverse shear force intensities. However, it is assumed that there 
are no transverse shear strains in this membrane theory and hence Q\ and Q 2 are 
not directly related to displacement variables (see Eq. (74b)). Hence, Q\ and Q 2 
need to be obtained from the first two equations of Eq. (79) as: 

Qi — T~ {M\x + Mq v — M 2 k 4 — M§k $) , 

1H ; ei < 85 “> 

Q 2 — z — ; {M2y 4- Mq x 4- M§k 4 4- M\k 5) . 

i + e 2 

The third equation of Eq. (79) is a statement of the balance of internal moments 
with respect to the £ axis, which is given by: 

Nq — N§ + M 2 /t6 2 — Mifcei “I" M§k\ — Afg/c 2 — 0 . (85b) 

Hence, only the three equations in Eq. (82) and the boundary conditions in Eq. (80) 
need to be solved in the nonlinear analysis of doubly-curved membranes. Equation 
(82) can be rewritten as: 


^ N\ x + N§2y 
< N&lx 4- N2y 
, Q lx 4- Q2y 




kiQi - k^Nei 4- /C62Q2 - k 4 N 2 j 


fpx] 


' mil 1 

kftiQx 4- k$Ni 4- k 2 Q2 -1- k 4 N&2 ] 

>+[T}< 

P2 

>=[T}< 

mv 

k\N\ — k 6l N 61 — Aj62^62 “ ^ 2^2 J 

1 ! 

Ip3 J 

! 1 

[ mw J 


(86) 


where iV 6 i = Nq/{1 + ei) and N 62 = iV 6 /( 1 4- e 2 ). 


2.3.1. Without bending stiffness 

If 76 is small and the membrane thickness h is so thin that the change of strain 
values with z is negligible, the strains in Eq. (63) can be approximated as: 

#22=e 2 , 2 Bi 2 = 7ei 4- 762 • (87) 

Substituting Eq. (87) into Eqs. (72) and (73) reveals that Mi, M 2 , and M 6 should 
disappear from the governing equations. Hence, Eq. (85b) becomes a null statement 
of N$ = iVe, and Eq. (85a) becomes 


Qi =Q 2 = 0. (88) 

Substituting Eq. (88) into Eq. (86) yields 


f Nl* + N 62y ' 


—k^Nsi — k 4 N 2 


(pi] 


r mil ' 

1 Neix + ^2y 

r 

k$Ni 4- k 4 N 62 

>+{T}< 

P2 

>=[T\< 

mi) 

[ 0 J 

1 1 

[ —k\Ni - keiNsi - ^62^62 - k 2 N 2 J 


[p3 J 

1 1 

[ mw 


( 89 ) 







I 
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These are the only three governing equations, and the boundary conditions in 
Eq. (80) reduce to: 

Along x = 0, X: 

5u = 0 or NiTii 4- Nq\T 2 i , 

Sv = 0 or 7 ViT 12 + , 

— 0 or N\Tis + ^61^23 • 

Along y = 0, Y: 

Su = 0 or iV 62 r n 4 N 2 T 21 , 

= 0 or A 7 g 2 Ti 2 4 N 2 T 22 1 (90) 

= 0 or N e2 Tiz 4 N 2 T 2 $ • 

We note that, if x and/or y are curvilinear, the magnitude of the differential 
area dxdy may change with x and/or y. This area effect is not accounted for in 
using integration by parts in Eqs. (73) and (76) because it can be done only if the 
explicit expressions of x and y are known. Next we show how to account for this 
area effect in using Eqs. (85a) and (86) to obtain the fully nonlinear theory of thin 
circular plates. 

2.3.2. Circular plate theory 

Figure 2 shows a general axisymmetric membrane described by an orthogonal 
curvilinear coordinate system xyz, where x denotes the meridian (also called the 



Fig. 2. The two coordinate systems used in describing the undeformed geometry of an axisym- 
metric membrane. 
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generator) and y denotes the circumference. For a circular plate, if the x is chosen 
to coincide with r, we have 


dx = dr , dy — rdO . 


(91a) 


Then, one can follow Eqs. (l)-(4) to show that all initial curvatures are zero except 
that 

k° 4 = - . (91b) 

r 

It is clear from Eq. (91a) that, if d6 and dr are fixed, the differential area dxdy 
changes with r. To account for this area effect caused by curvilinear coordinates, 
one just need to do the following changes for the spatial derivatives of a stress 
resultant: 




Hence, the fully nonlinear 
(85a), and (92) as: 


dNirdOdr 1 


dr rdOdr 
dNirdOdr 1 


= N lr + N x k 0 4 , 
= N ie k° 4 . 


(92) 


rdO rdOdr 

circular plate theory can be obtained from Eqs. (86), 


' Ni r + NX + N 62e k° 4 ' 

| 

kiQi — ksNei + k^Q2 — k 4 N2 ' 


( Pi 1 

N 61r + N 61 k° 4 +N 20 k o 4 

+ 1 

keiQi 7 * k§N 1 + ^2^2 + k 4 N§2 

> + [T}< 

P2 > 

w Qlr + Qlk® -f- Q20^4 j 

1 

, —kiNi — keiN^i — ^ 62^62 — ^2^2 J 


Us , 


= [T}{ 


mu 

mv 

mii) 




(93) 


where 

Qi = — (Mir + Afi *2 + M* ek° 4 - M 2 k 4 - M 6 k 5 ) , 

(94) 

Q 2 = “ (M20/S4 7 - MQ r T M$k 4 + M§k 4 + M\ ks) . 

1 4- e2 


The boundary conditions are the same as those in Eq. (80) except that dx and dy 
need to be replaced with dr and rdO , respectively. Similarly, dx and dy need to be 
replaced with dr and rdO in the nonlinear expressions of ei, e 2 , 761? 762 ? &i, fo? and 
k§ (= ^6i + k$ 2 ) shown in Eqs. (6), (7), (11a), (lib), and (18). 


2.3.3. Accounting for 76 

We point out here that the governing equations (83) and (84) are valid and can be 
derived using the vector approach only if 76 is negligibly small. If 76 is not small, 
substituting Eqs. (43) and (44) into Eq. (72) yields 

sn =L (Nidei + N 2^ e 2 + iVg i^ 7 ei + 7^62^762 M\5ki + M 2 ^k 2 + M§5k§)dxdy , 
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where 

N\ = Ni cos 76i 4 iV 6 sin 7 6 i , 

N 2 = N 2 cos 762 4 sin 7 6 2 , 
iV 6 i == (1 4- 6i)(JVe cos 7 ei - N x sin 7 6 i) , 

Ne 2 = (1 + e 2 )(N 6 cos 7 6 2 - A 2 sin 752 ) . 

Then, one can follow the same procedures to derive the governing equations. 
However, Eqs. (19), (20), (24a), and (24b) need to be used for 5e 1 , <5e 2 , 6761 , and 
< 5762 * The derivation is straightforward, but the governing equations obtained are 
complicated by 76 . However, if finite elements are used to discretize the system, 76 
does not complicate the formulation too much. 16 

Large scientific inflatable membranes are usually constructed from specially- 
shaped flat gores using a doubler material at the seams or other methods, 
which leads to structural models that defy closed-form solutions. Moreover, as 
materials for scientific membrane structures may exhibit nonlinear viscoelastic 
response, 17 the [Q] in Eq. (56) may become functions of strains and strain 
rates. Hence, the finite element method is more appropriate for analyzing real 
membrane structures. Furthermore, nonlinear finite element equations can be easily 
solved using an arc-length iteration method, and they can be used to deter- 
mine bifurcation points when wrinkling starts and to study bifurcated equilibrium 
paths. 


3. Axisymmetric Membranes Under Axisymmetric Loading 

Rotationafly sy mm etric membranes are popular in scientific applications because 
they are easy to manufacture and it is easy to predict and control their inflated 
shapes. Next we derive the fully nonlinear governing equations of axisymmetric 
membranes subjected to axisymmetric loading. 

From the undeformed geometry of the axisymmetric membrane shown in 
Fig. 2 and the coordinate systems and unit vectors shown in Fig. 1 we obtain 
that 


P = ai a -f r sin 0i b 4- r cos 6i c , 

dx = \Jda 2 4 dr 2 = ±day/l 4 r 2 , dy = rdO , 

dP ±1 . ^ 

Ji = -3- = 7 ===== (la + r a sm 0i b 4- r a cos 0\ c ) , 

dx y/lT rf 

• dP n- - /r 

J 2 = — = cos 6\ b - sm 0i c , 

ay ±1 

J3 = ji X j 2 = 7 ■ — (-r a ia + sm 6i b + cos 9 \ c ) , 

Vl + r2 


(95) 


where P is the position vector of point A, r is a function of a, r a = dr /da, and 
the ± is determined by the sign of da/dx. Substituting Eq. (95) into Eq. (3) we 
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obtain 


rJl + r 2 


Jc° — a h° — Jc° — h° — fl 

^4 — /r— 7> » ^61 — k 62 ~ “ u ■ 


rwl + r 2 


Since x and 6 are two curvilinear Lagrangian coordinates to be used, 
one needs to account for the area effect by making the following 
replacements: 

xr dNirdOdx 1 . n 

dx 7dMi = Nl ’ + N,k <- 

Li/JU i Ufl/UU^ //v/n v 


dNirdOdx 11 
ly rdO rdOdx r 16 ' 

The same replacements need to be done for all other spatial derivatives of stress 
resultants in the governing equations. 

For a specific membrane geometry, one can use Eq. (96a) to derive initial 
curvatures. For example, if it is a parabolic membrane with the apex at the origin 
of abc and the focus at a = a, then 


*o = L .0 _ 

1 (r 2 + 4a 2 ) 3 / 2 ’ 2 


ry/ r 2 + 4a 2 


If it is a spherical membrane with the center at the origin of abc and the radius is 
R, then 

R 2 = a 2 4- r 2 , r a = — , r aa — , 


jlO _ jlO 1 L.0 

^1 — «; 2 — — , ac 4 


/i? 2 - r 2 
rR 


Under axisymmetric loading, displacements only happen on the xz plane, and 
the convected axes f and fj are always perpendicular to each other. Hence, we have 

^61 — &62 — — 0 , 

h = h = T 2 iji + T 22 j 2 + T 23 j 3 = j 2 , ( 98a ) 

ii = ii = Tnji 4 - T12J2 4 - T13J3 = Tu ji 4 - ^I 3 j 3 1 
i3 = i x x i 2 = — T13J1 + Tnjs * 

Moreover, because ij = ii and i 2 = i 2 , 

761 = 762 = 76 = B12 = 0 , 

(98b) 

£11 = ei + z(fci - fcj) , #22 = e 2 + z(k 2 - k %) . 
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Equation (98b) indicates that there is no in-plane shear straining in axisymmetric 
membranes undergoing axisymmetric deformations, and hence x and y are the 
principle strain axes. Moreover, because of Eq. (98b), one can see from Eq. (73) 
that N§ should disappear from the governing equations. Substituting Eqs. (96a), 
(98a), and (98b) into Eqs. (14a), (10), (16), (17), and (70) yields 

Tn 0 T13 , 

1 + u x + wk\ _ w x - uki 


[T] 


0 

-Tn 


0 

Tii 


I’ll = 


1 + e, 


T 13 = 


1 + e. 


(99a) 



'0 

0 



' 0 


0 ' 

[Ki] = 

0 

0 

0 

, \k 2 ] = 

— A *4 

0 

- -k2 



0 

0 _ 


0 

fc 2 

0 . 


(99b) 


Pi = -p(l + ei)(l + e 2 )Ti 3 , p 2 = 0, p 3 = p(l + ei)(l + e 2 )Tii. 


(99c) 


3.1. Forward analysis without bending stiffness 

Under static axisymmetric loading and without bending stiffness, because of 
Eqs. (96a), (98a), (98b), and (99a)-(99c), the second equation of (89) becomes 
a null statement of 0 = 0 and the other two equations become 

N lx — k±N 2 = muTu -f mwT \ 3 , (100a) 

-k\N\ - IC 2 N 2 + p(l + ei)(l + e 2 ) = -muTu + mwTn . (100b) 

Accounting for the area effect by using Eq. (96b) in Eq. (100a) yields 

(101a) 

Post- multiplying Eq. (16) by [T] and using Eqs. (99a) and (99b), we obtain that 


= iV 2 &4 - N\k 4 + milTn + mwTi 3 . 
ax 


^=r 13 (fci-fc?), 

(101b) 

^ =Tn(fc?-fc 1 ). 

(101c) 

Moreover, it follows from Eq. (99a) that 


^ = (l + e 1 )T n -mfc?-l, 

(lOld) 

= (1 + 6i)Ti3 + uk\ . 

(lOle) 


Substituting Eqs. (96a), (98a), and (99a) into Eqs. (7) and (18) we obtain that 

e 2 = uk° 4 + wk%, k 2 = T n l4-Ti 3 kl, k 4 = T 13 k% + T n fc$ . (102a) 
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Because 76 = 0 and the bending stiffnesses are neglected, it follows from Eq. (74b) 
that 


ei = 


— Q12C2 
Q11 


^2 = $ 21*1 + 022^2 • 


Furthermore, we obtain from Eq. (100b) that 

p{ 1 + ei)(l 4- e 2 ) - k 2 N 2 4- muTiz - mwT n 






(102b) 


(102c) 


Substituting Eqs. (102a)-(102c) into Eqs. (101a)-(101e) yields five equations in 
terms of A/i, Tn, T13, u, and In other words, Eqs. (101a)-(101e) are the five 
governing differential equations. However, it is a fourth-order system because it 
follows from Eq. (90) that the boundary conditions are to specify: 


u or N\Tu 
w or A/i7i3 


(103) 


at x = 0 and x = L, where L is the arc length of the meridian. Because i 1 is a unit 
vector, the constraint equation needed is: 


I?1 +2i3 = l- (104) 

Equations (101a)-(101e) (or Eqs. (101a), (101c)-(101e) with T n = ±y/l — Tf s ) 
will be used to obtain numerically exact solutions by using a multiple shooting 
method. 18 

One can see from Eq. (96a) that, at r = 0, k £4, k 2 , and k$ may become 
singular and it is impossible to compute e x , N 2 , and fc x using Eqs. (102b) and 
(102c) because e2, k 2 , and k± cannot be computed using Eq. (102a). However, 
because of the axisymmetric geometry, we have 


Tn — 1 , T13 — 0 , — &2 > k 2 ki , 

k 4 = k%, e 1 ~e 2 , N\ = N 2 , N ix = 0 . 


at r = 0, where is obtained from Eq. (102a) by using Tn = 1 and 7\ 3 = 0. 

Note that k 2 and k± do not represent actual curvatures and they are functions used 
for exact description of £ and 7, as explained right after Eq. (18). Substituting 
Eq. (105a) into Eqs. (102b) and (102c) yields 


ei = -s 


Ni 


Q n + Q 12 


N 2 = (Q21 4- Q22)ei , = 


p(l 4- e x ) 2 — mw 
2iVi 


(105b) 


We note that Ni = N 2 in Eqs. (105a) and (105b) results in Qn = Q22- In 
other words, if an initially axisymmetric membrane is not manufactured to have 
Q11 = Q22? the membrane will not have an axisymmetric shape under axisymmetric 
loading. 

If a spherical membrane undergoes uniform expansion with respect to its center, 
we have 


u = Ti 3 = 0 , Tn = 1 . 


(106a) 
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Because of symmetry, we have 

Ni = N 2 , k\ = k 2 ? fe? — A:® , w x = iVi x = 0 . (106b) 

Moreover, substituting Eq. (106a) into Eqs. (99a) and (102a) yields 

ei = wk\ = e2 = ^ , (106c) 

It 

k 2 = ^2 , &4 = fc® - (106d) 

Equation (106d) is due to the fact that k 2 and £4 axe defined with respect to the 
undeformed length dx and are not the actual deformed curvatures, as explained 
right after Eq. (18). Equations (106a)-(106d) satisfy exactly the governing equa- 
tions (101a)-(101e) and make them to be 0 = 0. Substituting Eqs. (106b)-(106d) 
into Eqs. (102b) and (102c) yields the nonlinear function of p in terms of w as: 

2 (Q 11 + <3i2)u> + R 2 mw 

P (R+w) 2 ' 1 ’ 

If w/R is small and the pressure is static, one can replace R + w with R and then 
p is proportional to w. 


3 . 2 . Forward analysis with bending stiffness 


If bending stiffnesses are to be included in the analysis of axisymmetric membranes 
under axisymmetric loading, the strains shown in Eq. (63) need to be used, instead 
of Eq. (87). However, the surface analysis shown in Eqs. (95)-(99c) are still valid. 
Although Eq. (74b) show that and M$ may be nonzero due to elastic couplings, 
Eq. (73) shows that N§ and M& should not appear in the governing equations 
because fcg — 7ei = 762 = 0 (see Eqs. (98a) and (98b)). However, M x and M 2 
should be included. Moreover, substituting Eq. (98a) into Eq. (85a) yields Q 2 = 0 
and 

^- = (l + e 1 )Q 1 + M 2 k 4 - MifeS , (108a) 

ax 

where M\ x is replaced with M Xx + M\k\ to account for the area effect, as shown in 
Eq. (96b). The first and third equations of (86) become 


-• = k±N 2 — hQi — k®N\ + milTii 4 mwTis , (108b) 


= kiNi -b k 2 N 2 - p( 1 4 ei)(l 4- e2) - k%Qi - muTis 4 mwTn , (108c) 

where N Xx and Q Xx are replaced with N Xx + N x k% and Q\ x 4- Qik%, respectively, 
to account for the area effect. Moreover, Eqs. (101b)-(101e) are still valid and are 
repeated here as: 


dT n 


= r 13 (k ! - fc°) , 


dx 


(108d) 
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dT is 
d x 


= T 11 (k° 1 -k i), 


(108e) 


— = (1 + ei)Tn - wki 


1, 


(108f) 


— = (1 + ei)Tj3 + uki . (108g) 

Moreover, Eq. (102a) is also valid and is repeated here as: 

62 — uk 4 -|- wkQ , ^2 =: ^11^2 — ^13^4 i &4 — ^13^2 "l - -^11^4 * (109a) 

Because 76 = ke = k$ = 0, it follows from Eq. (74b) that 

Qnei 4 - Qu(k\ — ki) — N\ - Q1262 — Qi 5(^2 — fc§) i 


Q 4 l 6 l 4 - 044 (fcl “ ) — Mi — Q4262 “ Q 45 (fe — ^2) • 

Using Cramer’s rule we obtain that 

i [Mj - Q4262 ~ < 345(^2 ~ ^2)]Qll ~ [M ~ Ql262 ~ Qlbfa ~ 

1 1 QuQu-Qit 

_ [Ni - Ql 2^2 - Qlb {^2 - ^2)1^44 - [Ml - Q4262 - ^45(^2 - k 2 )}Ql 4 

QuQ 44 “ Q14 

(109b) 

Moreover, we obtain from Eq. (74b) that 

N2 = Q2161 4 - Q2262 -I- ^24(^1 — k^) 4 - <325(^2 — ^2) » 

M2 = Q5161 4 - Q5262 4 - Qs 4 (^i — &?) 4 - < 355 (fe — ^2) * 

Substituting Eqs. (109a)-(109c) into Eqs. (108a)-(108g) yields seven equations in 
terms of Mi, N\, Q 1, Tn, T13, u , and w. In other words, Eqs. (108a)~(108g) (or 
Eqs. (108a)~(108c), (108e)-(108g) with T n = ±y/l - T? s ) are the governing dif- 
ferential equations. We need the constraint equation (104) because the system is of 
sixth-order, as shown by the following boundary conditions: to specify 

u or ATiTn - Q1T13 


(109c) 


w or iSTiTis + QiTn (110) 

Tu(i.e.,^2) or Mi 

at x = 0 and x = L. 

Because of the axisymmetric geometry, the singularities at r = 0 can be dissolved 
because, at r — 0, 

T11 = 1 ? Ii 3 = 0 , k^k®, k\ = k 2 , fc 4 = ^4 , ei = e 2 , 

Ni = N2 , Mi = M2 , Mix — Nix — Qi — 0 . 


(111a) 
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Using Eq. (111a) in Eq. (109b) we obtain 

Nl(Qi4 + Qio) - MliQli + Q15) 


ei = 


(Q 11 + Ql2)(Qi4 + Q 45 ) — (Q 14 + Ql5)(Q*l + Q 42 ) 
MjXQn + Q 12 ) - Ni(Qu + Q 42 ) 


(111b) 


k 1 =ki+ 77 . 


(lllc) 


(Qll + Q12XQ44 + Q45) “ (Ql4 + Ql 5)(041 + Q42) 

Using Eq. (111a) in Eq. (109c) we obtain 

N2 = (Q21 4 - Q22X1 4 - {Q24 4 - Q25X&1 — &?) ) 

M2 = (Qsx 4 - Q52X1 4 - (Qs 4 + Qtt)(k\ - k ?) . 

If the material is isotropic, it follows from Eq. (98b) that the elastic energy II 
is given by: 

"-i/.oy fi)*** 

^-4-f 1 1- 

*■ X J 


2n rdx, 


IQ] 


(112) 

Here the first term is due to in-plane stretching and the second term is due to 
bending. 

3-3. Eigenvalue analysis 

To derive linear natural frequencies and mode shapes of vibration with respect to 
a deformed static equilibrium configuration we assume 

u = u4fi, w — w + w , (113) 

where u and w denote large static displacements, and u and w denote small dy- 
namic displacements. If the dynamic displacements are assumed to be harmonic at 
a natural frequency a;, we have 

u = —U) 2 U, W — —L J 2 U). (114) 

Substituting Eqs. (113) and (114) into Eqs. (101a)-(101e) and (102a)-(102c) and 
using Taylor’s expansions yields the following first-order expansions: 

dN x 


= iV 2 fc4 + iV 2 fc 4 - Nik% - muuj 2 Tii - mwu; 2 Ti3 , 


dx 


(115a) 


dT u 


= Ti s (*i - fc?) + f 13 *i , 


dx 


(115b) 
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^ = T n (*J - fci) - fnfci , 

(115c) 

— = (14 ei)Tii 4 eiTii — ivk® , 

(115d) 

— = (1 4 ei)Ti3 4 eiTi3 4 uk® , 

(115e) 


e 2 = uk 4 4 wk 2 , k 2 = Thfc!? — T13/C4 , £4 = T13&2 + Tufe® , 


(H5t) 

(116a) 


ei 


Ni - Q\2&2 


On 


i^2 = 021^1 4 022^2 , 


(116b) 


h 


p(e 1 4 e 2 4- eie2 4- e^) — k 2 N 2 ~ k 2 N 2 4 mwoj 2 Tu — k\Ni 


(116c) 


where Eq. (115f) is based on the fact that the natural frequency u 0 is the same 
for every point on the membrane. At r = 0, it follows from Eqs. (113), (114), and 
(105b) that 


ei = -~ 


N x 


N 2 = Nx, 


h = 


On 4 O 12 
p(e 1 4 eiei) 4* (mwJ 2 j 2) - k\Ni 
Wx 


(117) 


Substituting Eqs. (116a)-(116c) (or Eq. (117)) into Eqs. (115a)-(115f) 
yields six equations in terms of Nx, Tn , T13, w, and u. Because 
(7h 4- Tn) 2 4- (fi3 4- T13) 2 = 1 and T fx 4- T^ = 1, the linear expansion 
of Tn is Tn = — T13T13/T11. Hence, only Eqs. (115a), (115c)-(115f) need to 
be solved in the multiple shooting process. 


3*4. Inverse analysis 

In the design of regular structures (e.g., building and aircraft structures), the de- 
formed shape is usually not specified but is put as an inequality constraint in the 
design optimization process. Hence, the final design can be obtained by trying 
different undeformed shapes and dimensions to satisfy all equality and inequality 
constraints, which is not really an inverse problem. On the other hand, a precision 
membrane is usually required to have a precise inflated (or deformed) shape with a 
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Fig. 3. The generators of the undeformed and deformed shapes of an axisymmetric membrane. 

known inflation pressure, but the undeformed shape is unknown. Hence, the design 
of a precision membrane is a fully inverse problem. Next we show how to solve such 
a fully inverse problem. 

Figure 3 shows how the deformed shape of an axisymmetric membrane is re- 
lated to its unknown undeformed shape. Because of numerical difficulties in using 
coordinates a or c (= r) as the independent variable, we will use the curvilinear 
coordinate £ as the independent variable. Because 

(1 + ei)dx = (1 + e 2 )dy = <//?, (118) 

it follows from Figs. 2 and 3 and the definitions of and shown in 

Eq. (18) that 


n da . da . . da 

h = di' kl = ei = (1 + ei) ae’ 

(119) 

k% = shxa , k 2 = Sm& ( l + ea). 
c c 

(120) 


Moreover, because N\ represents the tension force per unit of dy and N 2 represents 
the tension force per unit of dx, we obtain from Eq. (118) that 

= Ni{ 1 + e 2 ) , N 2 = N 2 ( 1 + ei ) , (121) 

where N\ represents the tension force per unit of dr] and N 2 represents the tension 
force per unit of d£. From the free-body diagram of the bottom part of the deformed 
membrane cutting by a plane parallel to the c axis (see Fig. 3), because the internal 
pressure is balanced by the membrane force N \ , we obtain the following equilibrium 
equation: 


pc 2 tt = JVi27tc sin a . 


( 122 ) 
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On the other hand, Eq. (100b) is an equilibrium equation for a differential mem- 
brane element at any location. Substituting Eqs. (1 19)— (121) into Eq. (100b) and 
assuming static deformation (i.e., u = w = 0) yields 

da - sin a ~ 

P=^ N 1 + N 2 . (123) 

uc, c 

Because the pressure and the deformed shape are assumed to be known and hence 
p, d(£), and c(£) are known, Ni(£) and JV 2 (f) can be solved from Eqs. (122) and 
(123) as: 


1 2 sin £ ’ , . 

« / «* \ ( 124 ) 

N 2 = -r-7 (p- -£ Ml) ■ 
sin a \ at, J 

It follows from Eqs. (121) and (56) that 

iVi(l 4- e 2 ) = Qiiei + Q\ 2^2 » 

(125) 

N 2 {1 4 - ei) = Q 12^1 + Q 22^2 • 

Hence, ei(f) and e 2 (0 can be obtained from Eq. (125) to be 
c iVlQ22-iV2(0l2-iVl) 

61 Q11Q22 - (Qi 2 - Ni)(Qn - JVh) ’ (126) 

N2Q11 - Nl(Ql 2 - n 2 ) 

62 Q11Q22 - (Q12 - tfiXQu - #a) ' 

Because the circumferential strain e 2 at r = c should be constant due to symmetry, 
we have 


e 2 = - - 1 . 
c 


Hence, c(£) can be obtained to be 


1 +e 2 


Moreover, it follows from Eq. (118) that 

1 dx 

Hence, x(£) can be obtained by performing the following integration: 


It follows from Fig. 3 and Eq. (118) that 


da = ± V dx 2 - dc 2 = ±\ 1 - ( — ] dx = ±\ 1 — ( (1 4- ei) — 


dcy d£ 


d£j 1 + d 
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where the sign is determined by the sign of sin a . Hence, a(f ) can be obtained by 
perfor min g the following integration: 


a = 



(132) 


where a(0) — 0 is assumed. After a(f ), c(f ), and x(f ) are obtained from Eqs. (132), 
(128), and (130), one can easily obtain the undeformed geometry for manufacturing 
use. We note that it is more convenient to obtain the dc/d£ in Eq. (132) by numerical 
differentiation. 


4. Numerical Results 


Large static deformation analysis, linear eigenvalue analysis, and nonlinear inverse 
analysis of several different membranes are performed and discussed in this section. 
For all cases we consider axisymmetric membranes made of Kapton film that has 


E = 5.51581 GPa (800 kpsi) , j/ = 0.3, p = 1390 kg/m 3 , 
h = 0.0127 mm (0.0005 in) . 


(133) 


4.1. Forward analysis 

Inclusion of bending stiffness in membrane analysis results in difficulties in solution 
convergence in the multiple shooting process because M\ and Q\ have very small 
values. Figure 4(a) shows the undeformed (broken line) and deformed (solid and 
dotted lines) shapes of an axisymmetric parabolic membrane (see Eq. (97a)) having 
a/L — 20, 0 < x < L — 1.5 m and h = 15 mm, being subjected to an internal 
pressure of 500 Pa, and having a fixed rim at x = L. The boundary conditions for 
Eqs. (108a)-(108c), (108e)-(108g) are 

Qi = T 13 = u = 0 at x = 0 , Mi=u = w = 0atx = L. (134) 

Although Fig. 4(a) shows that the bending stiffness does not cause dramatic change 
of the deformed geometry, but Fig. 4(b) shows that N\ and N 2 are significantly 
changed by the bending stiffness. Moreover, Figs. 4(c) and 4(d) show that Mi, 
M 2 , and Qi are also significant. However the elastic energy ratio R e (= bending 
strain energy/total strain energy) is obtained using Eq. (112) to be 10.6% for this 
thick thickness case. If a/L = 20 , p — 500 Pa, and h = 15, 10, 5, 3 mm, R e is com- 
puted to be 10.60%, 5.08%, 1.30%, 0.40%, respectively. We note that R e decreases 
significantly when h decreases, and it becomes more difficult to obtain a convergent 
solution when R^ < 1.0%. If a/L =10 , h = 10 mm, andp = 500, 400, 300, 200, 100, 
10, 1 Pa, R^ is computed to be 2.69%, 2.79%, 2.91%, 3.06%, 3.25%, 3.47%, 3.50%, 
respectively. We note that R e increases slightly when p decreases. If h = 10 mm, 
p = 200 Pa, and a/L = 20, 15, 10, 5, 1, R e is computed to be 6.03%, 4.39%, 3.06%, 
2.13%, 1.41%, respectively. In other words, R e decreases when a decreases. Hence 
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Fig. 4. The deformation of an axisymmetric parabolic shell with a = 20L, h— 15 mm, and p = 
500 Pa: (a) deformed geometries with ( — ) and without (• • •) bending stiffness (NBS = no bending 
stiffness), (b) distributions of N\ and iV 2 , (c) distributions of M\ and M 2j and (d) distribution 
of Q i. 




Fig. 5. The deformation of a circular membrane having R = 1.5 m and a fixed outer rim and 
being under different pressures: (a) deformed geometries and (b) distributions of N\ and jV 2 . 
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the bending effeet is negligible especially when h and/or a are small. Because we 
consider only membranes made of the thin Kapton film shown in Eq. (133), bending 
stiffness will be neglected in all following cases and Eqs. (101a), (101c)-(101e) will 
be used in the forward analysis. 


4.1.1. Circular membrane 

For a circular membrane having a radius R = 1.5 m and a fixed outer rim, Fig. 5 
shows the deformed shapes and the distributions of N% and JV2 corresponding to 
different pressures. The corresponding boundary conditions are: 

T13 — u — 0 at x = 0 . a — w = 0 at x = R. (135) 

The multiple shooting solution is obtained using 21 shooting points. For p = 
0.322 Pa, tu(0) is obtained to be 18.6573 mm, which agrees well with the 




(c) 



(d) 


Fig. 6. The deformation of axisymmetric membranes whose undeformed shapes are parabolic: 
(a) deformed and undeformed geometries with p = 1 Pa and different a/L, (b) distributions of N x 
and N 2 with p — 1 Pa and different a/L, (c) deformed and undeformed geometries with p = 100 Pa 
and a/L = 0.5, and (d) distributions of N x and N 2 with p — 100 Pa and a/L = 0.5. 
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MSC/NASTRAN solution, 18.6589 mm, using 50 membrane elements and the 
solution, 18.6610 mm, using a different formulation and an integration method. 9 We 
note that MSC/NASTRAN uses Green strains and second Piola-Kirchhoff stresses. 

4.1,2. Paraboloidal membrane 

Figure 6(a) shows the undeformed (broken line) and deformed (solid line) shapes of 
axisymmetric membranes having parabolic undeformed shapes with different focus 
lengths (i.e., a), 0 < x < L = 1.5 m, and a fixed rim at x — L, and being subjected 
to an internal pressure of 1.0 Pa. The boundary conditions are given by: 

T13 = u = 0atx — 0, u = w — 0aitx = L. (136) 

We note that, under the same pressure, the deformation w(O) increases with 
a and hence the deformed geometry deviates more from a paraboloidal surface. 
Figure 6(b) shows the distributions of the internal force intensities Ni (along £) 
and N 2 (along rj). We note that N\ > N 2 , N\ = N 2 around the apex, and the area 
where N\ is close to N 2 increases when a decreases. However, Figs. 6(c) and 6(d) 
show that N 2 can be larger than Ni if a is small and/or p is large. 


4.1.3. Toroidal membrane 


We consider a toroidal membrane whose undeformed cross-sectional radius is R 
(= 3/ 7r m) and the center of the cross section is at c = Co = 10 m, as shown in 
Fig. 7(a). It follows from Fig. 7(a) that 

R 2 = a 2 + (r — Co) 2 , a = R sin , r = Co + R cos ^ 

a = ^ , 0<x<L = 3m. 

R 2 

We note that the meridian of the toroidal membrane starts from the point (a, c) = 
(0, Co + R). It follows from Eqs. (137) and (96a) that 



r ° = qT r r’ raa 


ft 2 k o = L 

(Co - r) 3 ’ 1 R’ 


k% = 


r - Co 

Rr 



(138) 


Figures 7(b) and 7(c) show the deformed geometry of the right cross section 
and the distributions of Ni and N 2 when p = 3000 Pa. The RMS of the deviation 
of the deformed cross section from a fitted circle 0.2% of i?, and this ratio increases 
with p. In other words, when p increases, the deformed cross section deviates more 
from a circle. 


4 . 2 . Eigenvalue analysis 

To understand how pressurization changes the natural frequencies of a membrane 
structure we consider a parabolic membrane having L = 1.5 m and a = 6L and 
perform eigenvalue analysis using the method shown in Sec. 3.3. Figure 8(a) shows 



Natural Praquanelaa (Hz) 
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the influence of the internal pressure p on the first five natural frequencies. It 
is apparent that natural frequencies increase about quadratically with p and the 
increment rate is high when p is small. Figure 8(b) shows the third linear mode 
shape when p = 60 Pa. As the derivation is for axisymrnetric deformation, only 
axisymmetric vibration modes can be obtained here. To obtain nonaxisymmetric 
modes one needs to solve the equations derived in Sec. 2.3 using the finite element 
method or other methods. 

In the formulations shown in Sec. 3.3 the change of p due to the small displace- 
ments u and w was neglected and it was assumed that p = p + p « p. However, 
it is necessary to include p in the calculation of low-order natural frequencies 
(especially the first-mode frequency) of some types of membranes. For example, 
if a pressurized spherical membrane undergoes the breathing- mode vibration (i.e., 
the first axisymmetric mode) and the gas inside the membrane is assumed to be 
an ideal gas, one can use p(R + tt)) 3 47r/3 = p(R + w + w) 3 47r/3 to derive the linear 
approximation of p as p = p( 1 — 3 w/(R + tD)). Hence, one can use Eq. (107) to 
derive the first axisymmetric-mode frequency uq to be 


UJl 


1 2(Qn + Q 12 ) + p(R + w) 
mR? 


(139) 


If p is not included, — y (2(Qn + Q 12 ) — 2 p(R + w))/(mR 2 ), which is erro- 
neous. However, p is negligible for high-frequency modes because of their wavy 
shapes (e.g., Fig. 8(b)). 


4.3. Inverse analysis 

4.3.1. Spherical membrane 

To verify the inverse analysis method shown in Sec. 3.4 we consider a deformed 
spherical membrane whose deformed radius is R (= 3/ir m) and whose material is 
shown in Eq. (133). It follows from Fig. 3 that 

a = R( 1 - cos a) , c = i? sin d, d=4, = 4 , 0 < £ < 3 m . (140) 

V R d£ R ” s - 

If p ~ 10 000 Pa, the undeformed geometry (broken line) obtained from the de- 
formed geometry (solid line) is shown in Fig. 9(a). The exact solution of w(Rir) can 
be obtained from Eq. (107) with w = 0 to be 

vR 2 

w(Rn) = a(Rir) — a(R7r) = — — . (141) 

Q 11 +Q 12 

The inverse solution of w(Rir) is 0.04553136 m, which is less than the exact solution 
(141) by 0.067%. Figure 9(a) shows also that the obtained undeformed geometry 
fits perfectly with a circle (dots). Moreover, Fig. 9(b) shows that N 2 = as it 
should be. Furthermore, we use the obtained undeformed geometry and the forward 
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Fig. 9. The inverse analysis of a known deformed spherical membrane with L — Rtt = 3m and 
p = 10000 Pa: (a) deformed ( — ), inverse undeformed (---), and circle-fitted undeformed (• ♦ -) 
geometries, and (b) distributions of Ni and JV 2 - 


analysis method to compute the deformed geometry, and the obtained deformed 
geometry overlaps with the initial deformed geometry. This verifies the proposed 
inverse analysis method. 


4.3.2. Paraboloidal membrane 

Because of their broad applications in co mmun ication, axisymmetric membranes 
having a deformed paraboloidal surface are important membranes. Next we consider 
an axisymmetric membrane whose deformed shape is parabolic and is given by: 

c 2 = 4aa. (142) 


From Fig. 3 we obtain 


<% - Vdc 2 + da 2 = ^1+ dc. 


(143) 


It follows from Eqs. (142) and (143) and integration that f and c are related as: 

£ = J\Z 4 + p + a log (c + a\J 4 + - a log(2a) . (144) 

Because f is the known spatial coordinate, c(£) can be obtained from Eq. (144) 
using a nonlinear root solver and then a(£) can be obtained from Eq. (142). More- 
over, it follows from Fig. 3 and Eqs. (142) and (143) that 


dc 


cos a = 


2a 


da 


da 


4 a 


--2 


d£ y/A a 2 T c 2 ’ 


sm a = — = 


y/Aa 2 T& ’ (4a 2 +C 2 ) 3 / 2 * 

(145) 


Figure 10(a) shows the deformed (solid line) and undeformed (broken line) 
shapes of an axisymmetric membrane having a parabolic deformed shape and a 




62 P. F. Pai & L. G. Young 



Fig. 10, The inverse analysis of a known deformed parabolic membrane with L = 1.5m, a=lOL, 

and p = 200 Pa: (a) deformed ( — ) and inverse undeformed ( ) geometries, and (b) distributions 

of N\ and N 2 . 


focus length a = 15 m, 0 < £ < L(= 1.5 m), and being subjected to an internal 
pressure of 200 Pa. Figure 10(b) shows the distributions of the internal force in- 
tensities Ni (along £) and N 2 (along 77 ). We note that N 2 > Ni and they become 
equal at x = 0. Substituting Eq. (145) into Eq. (124) yields 

y i = |V4a 2 + c 2 , N 2 = I (yiE? + c 2 + . (146) 

It is apparent that Ni < N 2 except that Ni = N 2 at c = 0 (i.e., the apex). However, 
Fig. 10(a) shows that the rim diameter at £ — L needs to be increased in order 
to have the deformed geometry is parabolic, but it results in N 2 > N\. On the 
other hand, Fig. 6(b) shows that, if the rim diameter at x = L is fixed, it results 
in N 2 < N\ but the deformed geometry is nonparabolic. Hence, one can reduce the 
increase of the rim diameter at £ = L in Fig. 10(a) to make N 2 close to N\ and 
keep them almost constant at every location. Of course, the deformed geometry 
will deviate a little from being a paraboloidal surface. 


4.3.3. Cylindrical membrane 

Figure 11(a) shows a known deformed circular cylindrical tube with two spher- 
ical ends subjected to an internal pressure p = 2000 Pa. If L\ — 2.5 m and 
R\ — R 2 = 2 m, Figs. 11(b) and (11c) show the obtained undeformed geome- 
try and the distributions of e\ and e 2 , respectively. The undeformed geometry has 
a sudden change at the junction of the circular tube and the spherical end because 
the spherical end has a smaller e 2 . If the deformed geometry and p are fixed, N\ 
and N 2 are also fixed. Hence, increase of e 2 requires decrease of Q t j, as shown in 
Eq. (125). In other words, the thickness h of the spherical ends needs to be re- 
duced because Q n, Q 22 , and Q\ 2 are proportional to h as shown in Eq. (56). From 
the geometry shown in Fig. 11(a) and Eq. (124) one can show that da/d £ = 0, 
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(c) 

Fig. 11. The inverse analysis of a known deformed circular cylindrical tube with two spherical 
ends, L\ = 2.5 m, R\ = R 2 = 2 m, and p = 2000 Pa: (a) the circular cylindrical membrane, 
(b) deformed and undeformed geometries, and (c) distributions of e\ and e 2 . 


c/ sin d = J?i, and N 2 = 2Ni = pi?! for the cylindrical tube; and dd/d£ = I/R 2 , 
c/ sin d = R 2 , and iV 2 — iVj = pR 2 /2 for the spherical ends. Replacing N 2 and 
Ni with pR 2 /2, Qij with qQ t j, and e 2 with the circumferential strain e 2c of the 
cylindrical tube in Eq. (126) yields 

J2//S r\ , q(Qi2 ~ Q11 + 2Qi 2 e 2c )pi? 2 (1 + e 2c )p 2 il| 

q {Q11Q22 - Q 12 )e 2c + ^ = 0 » 

(147) 

which can be solved for the q (= h n ew/horigmai) required to match e 2 at the junction. 
Figure 12(a) shows that, if the thickness of the spherical ends is reduced to be 
hnew = 0.005451 mm (i.e., q = 0.4292 from Eq. (147)), the undeformed geometry 
is smooth at the junction. 

Replacing N 2 and N\ with pR 2 / 2 and e 2 with the circumferential strain e 2c of 
the cylindrical tube in Eq. (126) yields 

+ «,(<?»- <3,2 + (<a _ Q n Q^ _ 0 , 

(148) 
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c(m) c(m) 

(a) (b) 

Fig. 12. The inverse analysis of a known deformed circular cylindrical tube with two spherical 
ends, L\ = 2.5 m, R\ =2 m, and p = 2000 Pa: (a) deformed and undeformed geometries if 
/?2 —2m and the thickness of the spherical ends is reduced to 0.005451 mm, and (b) deformed 
and undeformed geometries if the radius of the spherical ends is increased to i ?2 = 4.6597 m. 

which can be solved for the R 2 required to match e <2 at the junction without chang- 
ing the thickness of the spherical ends. Figure 12(b) shows that, if R 2 is increased to 
4.6597 m, the undeformed geometry does not have a sudden change at the junction. 

5. Concluding Remarks 

A fully nonlinear membrane theory is derived by treating membranes as shells and 
using Jaumann strains and stresses and new concepts of local relative displacements 
and orthogonal virtual rotations. Fully nonlinear governing equations for different 
membranes are listed, especially the ones for axisymmetric membranes subjected to 
axisymmetric loading. The membrane theory is used to obtain numerically exact de- 
formed shapes of several axisymmetric membranes subjected to internal pressures. 
Numerical results show also that bending stiffness is negligible for thin membranes, 
and natural frequencies of pressurized membranes increase with the internal pres- 
sure. Moreover, an inverse solution method for computing the undeformed geometry 
of a membrane with a known/designed deformed geometry and a known inflation 
pressure is presented and numerical results are obtained and verified. 
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